import cmath
a,b,c = 1,2,20
delta = b*b-4*a*c
x1 = (-b+cmath.sqrt(delta))/(2*a)
x2 = (-b-cmath.sqrt(delta))/(2*a)
print(x1,x2)

import sympy
x = sympy.Symbol('x')
rslt = sympy.solve(x**4+2*x+20)
print(rslt)

